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O ■ Abstract 

o : 

. We investigate classical formation of a D-dimensional black hole in a high energy collision of two 

CIhI particles. The existence of an apparent horizon is related to the solution of an unusual boundary- 
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value problem for Poisson's equation in flat space. For sufficiently small impact parameter, we 
construct solutions giving such apparent horizons in D = 4. These supply improved estimates of 
the classical cross-section for black hole production, and of the mass of the resulting black holes. 
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' We also argue that a horizon can be found in a region of weak curvature, suggesting that these 

o ■ 

solutions are valid starting points for a semiclassical analysis of quantum black hole formation. 
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I. INTRODUCTION 



The proposal that the fundamental Planck mass could be as low as a TeV has excited 
new interest in the problem of black hole formation in ultra-relativistic collisions. TeV- 
scale gravity scenarios offer a completely new perspective on the hierarchy problem, and 
arise via either large extra dimensions |jl|, or compact dimensions with large warp factor; 
a model for the latter appears in and string solutions in [Q. It has long been believed 
that high energy collisions where the center of mass energy substantially exceeds the Planck 
mass would produce black holes; this statement can be thought of as a simple extrapolation 
of Thome's hoop conjecture[^] and has been more recently discussed in [^. Lowering the 
Planck scale to OiTeV) thus raises the exciting prospect that black holes can be produced 
at future accelerators, perhaps even at the LHC|0, H-^ 

Clearly we would like to better understand this process. One important problem is to 
estimate the cross-section for black hole production. It has been argued that at high energies 
black hole production has a good semiclassical description P, 0], since in such cases a horizon 
should form in a region where the curvature is weak and quantum gravity effects are small. 
This leads to the naive estimate that the cross-section for black hole production is roughly 
given by 

a ~ ^^rl{^/s) (1) 

where denotes the horizon radius corresponding to CM energy -^i. One would like to 
make this estimate more precise to improve experimental predictions, and in particular to 
derive a differential cross-section depending on the mass and spin of the black hole produced. 
Furthermore, validity of the estimate (|I]) has been challenged |[10|, |Tlj|. In particular, while 



Penrose [ p!^] and D'Eath and Payne|13, 14, [T^ (for a comprehensive treatment see |16]) have 



studied the problem of a classical collision with zero impact parameter and shown that a 



closed trapped surface forms, |jTT[ argues that such collisions will not form black holes at 
non-zero impact parameter. Finally, collisions of cosmic rays with our atmosphere have 
energy reach beyond that of LHC, and reasonable conjectures about neutrino fluxes suggest 
the possibility that their black hole products might be seen with present or future cosmic 
ray observatories O, ITHl O, POjP or absence of their observation could place improved 



^ For a review, see 
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bounds on the fundamental Planck scale pT|. However, unlike the collider setting, where 



production above threshhold is copious, conclusive statements here depend sensitively on 
the exact factors in ([1|). This calls for improved estimates. 

This paper will revisit the classical problem of black hole production in high energy 
collisions. In particular, understanding the case of non-zero impact parameter is crucial 
to improving the estimate (|T]). We will use the methods of Penrose |T2| and D'Eath and 



Payne ||13|, |14|, |T^, where each incoming particle is modelled as a point particle accompanied 
by a plane-fronted gravitational shock wave, this wave being the Lorentz-contracted longitu- 
dinal gravitational field of the particle. At the instant of collision the two shock waves pass 



through one another, and nonlinearly interact by shearing and focusing. Penrose] 12] and 



D'Eath and Payne p3|, |T^, |T5[ studied the case of zero impact parameter b, and by finding a 



closed trapped surface, derived a rigorous lower bound of M > a/s/2 and improved estimate 
M ^ .84y^ for the mass of the resulting black hole. This paper extends this analysis to 
b ^ 0. In the next section we review the basic approximations involved, and generalize the 
construction of the incoming waves to arbitrary dimensions. We then reduce the problem of 
finding a closed-trapped surface to a rather unusual boundary value problem for solutions 
of Poisson's equation; this problem has a close analog in the physics of soap bubbles. Us- 
ing conformal methods which apply only in four dimensions, section three gives an explicit 
solution to this problem in that case and gives a lower bound 

C"bh production > 32.5(G'^s/4). (2) 

on the classical black hole cross section, where G is Newton's constant. Section four dis- 
cusses issues in extending this analysis to higher dimensions, and in better justifying the 
semiclassical approximation; we follow with conclusions. 



II. TRAPPED SURFACES IN SHOCK- WAVE GEOMETRIES 

In the TeV gravity scenarios in question, we will assume that the particles of the standard 
model propagate on a brane, but that gravity propagates in a higher- dimensional space 
with either large volume or large warping.^ This lowers the fundamental Planck scale Mp, 

^ For a brief unified review of these scenarios, see M. 
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perhaps to 0{TeV), while maintaining the observed value of the four-dimensional Planck 
mass, M4 ~ lO^^GeV. We expect to be able to create black holes in scattering processes with 
center of mass energy E > Mp. In general such black hole solutions will have complicated 
dependence on both the gravitational field of the brane, and on the geometry of the extra 
dimensions. However, there are two useful approximations [0] that may be used for a wide 
class of solutions. First, the brane is expected to have tension given by roughly the Planck 
or string scale, and so for black holes substantially heavier than the Planck scale the brane's 
field should be a negligible effect. Secondly, if the geometrical scales of the extra dimensions 
(radii, curvature radii, variational scale of the warp factor) are all large as compared to 
1/Mp, then there is a large regime where the geometry of the extra dimensions plays no 
essential role. Therefore it is often a good approximation to consider high-energy collisions 
in D-dimensional flat space. 

In the center of mass frame, each of the high-energy particles has energy 

/i = . (3) 

We use a coordinate system {u, v, x*) where retarded and advanced times {u, v) are {i—z, i+z) 
in terms of Minkowski coordinates, z being the direction of motion, and x\ i = 1 . . . D — 2, 
denotes transverse coordinates. The impact parameter is b and the particles are initially 
incoming along (a;*) = (±6/2, 0, . . . ,0). 

The gravitational solution for each of the incoming particles can be found by boosting the 
rest-frame solution. The gravitational field outside a particle is well approximated by the 
Schwarzschild solution. Since we will be concerned with long-distance phenomena such as 
formation of large horizons, short-range modifications of the solution should not be relevant. 
The D-dimensional Schwarzschild solution with mass M is 

where (if2|,_2 are the line element and volume of the unit D — 2 sphere. The 



Aichelburg-Sexl solution is found by boosting this, taking the limit of large boost and 
small mass, with fixed total energy ^. The result for a particle moving in the +z direction 
is the metric 

ds'^ = -dudv + dx''^ + ^{x')6{u)du'^ . (5) 
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Here $ depends only on the transverse radius p = \^x^Xi, and takes the form 



IGiiGp 

Note that $ satisfies Poisson's equation 

= -IGTcGpS^-^: 



D>A 



(6) 
(7) 



(8) 



in the transverse dimensions (here V is the D — 2-dimensional fiat-space derivative in the 
(x*)). This spacetime solution is manifestly fiat except in the null plane m = of the shock 
wave. If we consider an identical shock wave travelling along -u = in the —z direction, by 
causality these will not be able to infiuence each other until the shocks collide. This means 
that we can superpose the two solutions of the form (H) to give the exact geometry outside 
the future light cone of the collision of the shocks. 

The coordinates u, v, suffer the drawback that geodesies and their tangent vectors 
appear discontinuous across the shock. This can be remedied by going to a new coordinate 
system defined by 



u = u 



V = V + M{u) + 
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= x' + -VMx)e{u) 



(9) 
(10) 

(11) 



(where 6 is the Heaviside step function) in which both geodesies and their tangents are 
continuous across the shock at m = 0. In these coordinates, the metric of the combined 
shock waves becomes 



where 



H\f =5., + ^V,V,$(x-xi)m^(m) 



-du dv + 



rr(l) rr(l) , 7t(2) 7t(2) e 
^ik ^jk + ^ik ^jk - ^ij 



H. 



(2) 
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(12) 

(13) 
(14) 



with $ given by eqns. (|D and (^, and with 



xi = (+6/2,0,...,0) , X2 = (-6/2,0,... ,0) . 



(15) 



Here x = (x*) in the transverse flat D — 2-space. 

A marginally trapped surface 5 is a closed spacelike D — 2-surface, the outer null normals 



of which have zero convergence ||23|| . For the case of D=4 and 6 = 0, Penrose |p!2| found such 
a surface in the union of the two shock waves. This consisted of two flat disks with radii 
Pc dX i = —AGfilnpc, z = ±4G/ilnpc- Matching their normals across the boundary, which 
lies in the collision surface u = v = 0, then determined pc = ^Gp = r^- This construction 
immediately generalizes to the case D > 4, 6 = where 

= [a^J ■ '''' 

Generalizing to 6 7^ and arbitrary dimensions, we attempt to construct S in the union 
of the incoming null hypersurfaces v <{] = u and u < = f . These hypersurfaces intersect 
each other in the D — 2-dimensional surface m = = f , and S will intersect this D — 2-surface 
in a closed D — 3-dimensional surface C, to be determined. In the flrst incoming null surface 
f < = M, 5 will be defined by 

?; = -^i(x) with ^1 > interior to C, ^1 = on C, (17) 

and one may straightforwardly check that the outer null normals will have zero convergence 
for t> < as long as 

V^(*i - $1) = interior to C. (18) 
Similarly, in the second incoming null surface u < = 5 will be defined by 

u = -^2(x) with ^2 > interior to C, ^2 = on C, (19) 
and the outer null normals will have zero convergence for m < as long as 

V^(*2 - $2) = interior to C. (20) 

Finally, the outer null normal to S must be continuous at the intersection m = = f ; if not 
there would be a (5-function in the convergence. A necessary condition for continuity is 

V^i ■ V^2 = 4 on C ; (21) 

since \&q,, a = 1, 2, vanish on C, Vj^I/q, is normal to C and this condition is also sufficient. 
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Finding a marginally-trapped surface therefore reduces to a simple mathematical prob- 
lem. Specifically, note that (||), (|T8|) imply that \1/q, satisfies Poisson's equation with sources 
at x„. Define the rescaled functions 

satisfying 

V^^7(x, X,; C) = -Qd-^6''-\^ - x«), (23) 
^(x, x„; C) =0 for X on C . (24) 

These are thus the Dirichlet Green's functions for sources at xi,X2 and with boundary C. 
The problem of finding the marginally trapped surface is equivalent to the following: 

Problem C. Given two points Xi and X2 in Euclidean D — 2-space, and a constant B > 0. 
Let g{yi, Xq,; C) be the Dirichlet Green's functions satisfying Find a closed D — 3- 

surface C enclosing the points with the following property: 

Vx5(x,Xi;C) ■ Vx(7(x,X2;C) = 52 (25) 

for all points x on C 

As a trivial example, if xi = X2, then the unit D — 3-sphere about xi is a solution C to 
Problem C for i? = 1. This reproduces Penrose's trapped surface in D = A with suitable 
scaling, and gives its generalization, Eq. ([T6|) , for D > 4. Given general xi and X2, does a 
solution C always exist? Clearly not if the points are too distant from each other, because a 
collision at large enough impact parameter cannot produce a black hole. Given general xi 
and X2, is the solution C unique? We shall see that it is not. We also remark that solutions 
for different B are related by simple scale transformations. 

One way to understand Problem C is via another physical problem that serves as a 
simple analog. Consider a ring of wire with shape C in the x, y plane in three dimensions, 
and suppose that this ring is spanned by a soap film. If we apply pressure to the soap film, 
then in the limit of small displacement, its vertical displacement z{x, y) satisfies the equation 

V z{x,y) = , (26) 

a 

where a is the film's tension. Generalize to the case of two films, held apart by pressures in 



the +z and —z directions. If the pressure is exerted at points x^,, then the solutions to (|26D 



are the above Dirichlet Green's functions. If the horizontal positions of the pressure points 
are the same, and the ring C is circular, then the angles 6i, 62 of the soap films with the x, y 
plane are constant around C. Now separate the pressure points slightly in the x-direction - 
this will change these angles, and they will be functions of position along C. Problem C is 
that of finding the curve C for which 

52 = tan 01 tan ^2 (27) 

is constant over C, and equal to the value for zero x-separation of the pressure points. 
This problem can also be generalized to a higher-dimensional analog. One can argue for the 
plausibility of a solution, at least for small enough x-displacement, by noting that deforming 
a point on C towards the center increases the angles and thus the local value of and 
deforming away decreases B. This suggests that C can be adjusted point-by-point to make 
B equal to the given constant value over the ring. The next section will solve this problem 
explicitly in the special case D = A. 

III. EXPLICIT CONSTRUCTION FOR D = A 

In D = 4, Problem C may be readily solved, at least for sufficiently close points, by a 
trial- and-correction construction, in two steps: First, choose trial points Xi and X2, and a 
trial curve T; then construct 5f(x, x'; T), and evaluate 

/(5;T) = Vx(7(i,ii;T)- Vx(7(i,i2;T), 5 on T. (28) 

Second, correct the trial solution by finding a conformal transformation x = x(x) that sends 
/(x; T) to B^] this will send Xi and X2 to some points Xi and X2, and will send T to a curve 
C obeying eq. (pSf) . Thus we obtain a solution of Problem C. This works because Poisson's 
equation in dimension 2 is conformally invariant, while / transforms like 

'/(i;r). (29) 

Let us now construct some solutions to Problem C. Given two trial points Xi and X2 
such that I |xi — X2I = a < 1; we may take 

ii = (a,0), 52 = (-a,0). (30) 



/(x;C) 



9x 
9x 
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For the trial curve T choose the unit circle |x| = 1. Then the Green function evaluated for 
the points Xi and X2 is 

,(x,x,;r) =-2l-(^^-3Tj^^j> (31) 
,(x,x.;r) --^^-[ ^,^^^y^,.y. )^ (32) 

and 

/(x; T) = ^ . (33) 

To conformally send /(x;T) to S^, we use complex analytic variables z = x + iy and 
z = X + iy. We evidently require a transformation z{z), analytic on and within the unit 
circle T parametrized by 5 = exp(i0), with 

dz 



dz 



/(x; r) (1 _ 2ax + a2) (1 + 2ax + a?) ^^^^ 



[I - a') 



2\2 



(1 - a2exp(2i0))(l - a2exp(-2i(/))) 

By inspection, the required tranformation satisfies 

dz 1 — 

di ^ B{1 - a^P) 

which integrates to 



(35) 



(36) 



1 — , f 1 + az\ 

z(z) = In . 37 

^ ^ 2Ba \l-az ' ^ ^ 



The true points (xi,X2) are then found to lie at 



25a V 1 - a 

Restoring physical dimensions, 



1 - / 1 + a" , , , 

^ = ±^77^1n ^ o],y-0- (38) 



^1 =8G'/x^(x,xi;C), (39) 



*2 =8G'/iy(x,X2;C), (40) 
1 

4G^' 



B =:^, (41) 



f2Gii{l- a'^) ^ /l + a2 



6(a) = inl^] . (43) 



Thus we have constructed a marginally trapped surface S for any value of impact parameter 
b{a) that can be obtained from the above formula for some a, with < a < 1. The area of 
S is found to be 

Area(>S) = WniGfi)' ^^ 'f^' In ([^) • (44) 

Now S may or may not be an apparent horizon, which is defined as the outermost marginally 
trapped surface. However, the existence of S means either that S is in fact the apparent 
horizon, or an apparent horizon exists in the exterior of S. Because S can be shown to 
be convex, and because the 2-metric is Euclidean, Area(5) is a lower bound on the area of 
the apparent horizon. Modulo technical issues about cosmic censorship, the existence of an 
apparent horizon means that the collision will produce a black hole (or more than one black 
hole, though this seems unlikely in the present setup). Moreover, by the Area Theorem, the 
mass of the final black hole (assumed single) is bounded below. 



Mfinaibh > Jin ( YZr^ ) ' 



and the fraction of total energy 2/i = ^/s emitted as gravitational radiation is bounded 
above, 

^::^<i-i^jinfi±4y (46) 

In fact, -Egravrad niay be significantly smaller because the final black hole is expected to be 
rotating (unless 6 = 0), tying up substantial energy. 

The function h{a) (for given /i) reaches a maximum value of 

6„,ax ^3.219G/i (47) 

at 

a^,, ^ 0.6153 (48) 

so this is the greatest impact parameter for which this construction can demonstrate pro- 
duction of a black hole. The corresponding lower limit on the cross-section is 

abhproduction > Trfo^x ~ 32.552(G'/x)2. (49) 

Previous estimates for the black hole production cross-section used (|l|), which gives a ~ 
50.27(G/x)^, where Vh = 2Gy/s = AGfi is the Schwarzschild radius belonging to the total 
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energy available. Our lower limit is about 65% of this estimate. Another interesting quantity 
is the mass of the final black hole. Equation ( ^Sf ) together with ( ^3D gives a lower bound on 
the mass as a function of impact parameter. We find a range from 0.71^/s for 6 = to OAby/s 
for b = &max- The perturbative analysis of |T3|, [1^, ^ raises the former to M ^ .84y^; we 



expect a corresponding increase in the latter upon further analysis. 

At values a > Omax (but still a < 1) our construction produces a second, smaller 
marginally trapped surface for the same b. This shows that solutions to Problem C are 
usually not unique. 

It is also interesting to better understand the shape of the curve C for the D=4 solution. 
This is readily found from (|37D, and takes the form 



(1 - si^h^ r + (1 + sin^ r = a' . (50) 
4:Gfi{l - a^) 4G/i(l - a^) ^ ^ 

Obviously this approaches a circle of radius 4(7/1 as a — 0. Fig. |l] displays the curve C in 
the transverse collision plane for various values of b. 



IV. FURTHER RESULTS AND DIRECTIONS 

Clearly it would be desireable to find an explicit solution to Problem C in the higher- 
dimensional case, in order to give more careful estimates of the cross-section in the physically 
interesting situation where the extra dimensions are relevant. Nonetheless, we have given a 
heuristic argument for the existence of such a solution, which is buttressed by the explicit 
construction of the trapped surface in the four-dimensional case. This appears to demon- 
strate classical black hole formation at non-zero impact parameter — answering the criticism 



of 111 



For head-on collisions, 6 = 0, in D > 4, the apparent horizon is the union of two balls of 
radius Pc = fh, Eq. (|16]). The corresponding lower limit on black hole mass is displayed in 
Table |. 

Solving Problem C in D > 4 for 6 > may require numerical work. Note that by 
symmetry it still remains a two-dimensional problem of finding a curve C which produces 
the surface C as a surface of revolution about the axis connecting the two source points. 
However, in the D > 4 case, the relevant Green function is of the form (|^ and no longer 
transforms nicely under conformal transformations. Of course even solving this problem is 
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FIG. 1: The intersection curve C of the marginahy trapped surface S with the colhsion plane 
(n = = u). Several curves C, for various impact parameters b, are superposed; spacetime 
dimension is D = 4. Distances are in units of G^/s = 2Gfj,; in these units = 2. Incoming particle 
pairs appear in the horizontal line = at pair separation b; wider pairs therefore correspond to 
smaller curves C. Values of b are 0+, 0.6 x, 1.2*, 1.40, 1.55», and 1.609o, the last being maximal. 

only a starting point - it provides a lower bound on the mass of the black hole, but as was 



found in the axial symmetric case in |1^, the resulting black hole will absorb more 



of the energy contained in its external fields, and it should be possible to raise this bound 
by studying this subsequent evolution. Recall that in the axial-symmetric case in D = 4 
this resulted in an additional enhancement of approximately 119%. Going beyond this, the 
ultimate goal of such an analysis - and its quantum extension - is to compute the differential 
cross-section depending on the mass and spin of the resulting black hole. 
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TABLE I: Lower limit on M^h in head-on collisions 6 = 0, for various dimensions D. 
Here ^sch is the {D — 2-dimensional) horizon area of a Schwarzschild black hole of mass ^/s^ 





while j4trap is the area of the Penrose 


marginally trapped surface. 




D 


^trap /^Sch 


MBH/^/^ 


1 - Mbh/V^ 


4 


0.50000 


0.70711 


0.29289 


5 


0.54270 


0.66533 


0.33467 


6 


0.55032 


0.63894 


0.36106 


7 


0.55080 


0.62057 


0.37943 


8 


0.54928 


0.60696 


0.39304 


9 


0.54720 


0.59642 


0.40358 


10 


0.54502 


0.58798 


0.41202 


11 


0.54293 


0.58105 


0.41895 



One might also wonder if a black hole is produced in the high energy collision of a 
particle with a purely gravitational shock wave, in any dimension D > A. The present 
analysis suggests that the answer is "no". If we attempt to replace $2(x) by a source- free 
solution of Laplace's equation, to model a purely gravitational shock wave, then Eqs. (p!9| , pOD 
have no solution at all, by the maximum principle for elliptic equations. Therefore no 
apparent horizon exists in the incoming wavefront surface. Similarly, the collision of two 
purely gravitational shock waves seems not to produce a black hole. Moreover, in the 
collision of two particles at 6 > 0, as studied here, one might wonder if additional, smaller 
marginally trapped surfaces might appear, enclosing one particle but not the other; a similar 
argument shows not. We have not considered trapped surfaces that might exist to the future 
of the incoming wave surface, however, so these arguments cannot conclusively rule out black 
holes. 



It is also interesting to compare the estimate (^91) to a heuristic argument presented in 21 



Anchordoqui et. al. argue that one may improve estimates of the cross-section by taking 
into account angular momentum dependence of the Schwarzschild radius. Specifically, for 
CM energy ^/s and impact parameter b, the angular momentum is J = h^jl. One expects 
that the maximum impact parameter occurs near a value of h that equals the corresponding 
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angular-momentum dependent radius r/j. This is given byp5 



If we set 6 = in the former expression, that gives 



„_3 _ IQtxGM 
' h 



This leads to a cross-section estimate 

.2_ r. , iD-2r 

h 



{D - 2)^ 
16 



(52) 



2 

D-3 



'''"'^h.spherical • (^3) 



16 

In D=4, this gives a relative correction factor of 64%, closely matching the factor in 
The correspondence of these results suggests that this estimation technique may indeed be 
approximately correct in higher dimensions. 

A final point concerns the validity of the semiclassical approximation. We expect that the 
semiclassical approximation for black hole formation should be justified if a horizon forms 
at small curvature.^ Of course the surface S that we have constructed lies in the plane of 
the shocks, and thus is in a sense close to a region of large curvature; understanding the 
importance of this would require regulating the solution taking into account finite size/mass 
effects. However, we believe that this is not necessary, as it is possible to find a trapped 
surface outside of the planes of the incoming shocks. First, in the case 6 = 0, consider 
Penrose's flat disk S of radius p = pc in the incoming null wavefront surface m = 0. Construct 
the null plane TV emanating from S in the opposite direction v = —const; Af has zero 
convergence because it is a null plane. Deform S into the future along Af, while leaving it 
fixed in some neighborhood of its boundary C. This deformed D — 2-surface S' will still have 
zero convergence along the f-direction everywhere, and thus will be an apparent horizon, 
now with weak spacetime curvature on it. (In fact, S' has exactly the same area as S, and 
so gives exactly the same bound on black-hole mass.) 

Second, in the case 6 > 0, we can proceed similary. We can still construct a null surface 
TV emanating to the future from our original marginally trapped surface S, generated by null 
geodesies normal to S. However, TV will not be a fiat null plane, because the null normals to 



^ In the context of string theory, the Schwarzschild radius must also be larger than the string length 
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S have nonzero shear. We can still deform S some distance to the future along TV to create 
a new D — 2-surface S' but the convergence of the null normals of S' will go positive, due 
to the shear (we cannot go too far or we will run into a caustic, i.e. the convergence will go 
to +00). Thus S' will actually be a trapped surface, not a marginally trapped surface. This 
would mean that a marginally trapped surface must lie somewhere outside of it. In any case 
S' has weak spacetime curvature everywhere on it, and implies the existence of a black hole. 

V. CONCLUSIONS 

The existence of a closed trapped surface in the collsion geometry of two ultra-relativistic 
particles clearly demonstrates classical black hole formation. The argument that these sur- 
faces are present in the weak curvature region further suggests that this process can be con- 
sistently treated in a semiclassical analysis, and should help lay the foundation for a more 
rigorous justification of such an analysis. Furthermore, we have found improved estimates 
on the production cross-section for black holes. While these estimates are not enormously 
different from the more naive estimates of |]^, |^ , it is important to know their size in improv- 
ing discussion of the sensitivity of cosmic ray observations to black hole production. Future 
directions include more complete analysis of the higher- dimensional classical problem, which 
may require numerical work. The present classical analysis should serve as a starting point 
for a more complete investigation of the semiclassical approximation to black hole formation. 

After completing this work we were kindly informed by R. Penrose that he had found 
related unpublished results on the existence of a maximal impact parameter in D = 4, many 
years ago. 
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